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ABSTRACT 

emission is the dominant cooling mechanism in Jupiter's thermosphere 
and a useful probe of temperature and ion densities. The Hjj' ion is predicted 
to form in the thermospheres of close-in 'hot Jupiters' where its emission would 
be a significant factor in the thermal energy budget, affecting temperature and 
the rate of hydrogen escape from the exosphere. Hot Jupiters are predicted to 
have up to 10^ times Jupiter's emission because they experience extreme 
stellar irradiation and enhanced interactions may occur between the planetary 
magnetosphere and the stellar wind. Direct (but unresolved) detection of an 
extrasolar planet, or the establishment of useful upper limits, may be possible 



^Visiting Astronomer at the Infrared Telescope Facility, which is operated by the University of Hawaii 
under Cooperative Agreement No. NCC 5-538 with the National Aeronautics and Space Administration, 
Office of Space Science, Planetary Astronomy Program. 

^NRC Research Associate 

•^Partially supported by the National Aeronautics and Space Administration through the NASA Astrobi- 
ology Institute under Cooperative Agreement No. NNA04CC08A issued through the Office of Space Science. 



- 2 - 



because a small but significant fraction of the total energy received by the planet 
is re-radiated in a few narrow lines of within which the flux from the star is 
limited. 

We present the observing strategy and results of our search for emission from 
the Q(1,0) transition of Hjj' (3.953 fj,m) from extrasolar planets orbiting six late- 
type dwarfs using CSHELL, the high-resolution echelle spectrograph on NASA's 
Infrared Telescope Facihty (IRTF). We exploited the time-dependent Doppler 
shift of the planet, which can be as large as 150 km s"^, by differencing spectra 
between nights, thereby removing the stellar photospheric signal and telluric 
lines. We set hmits on the emission from each of these systems and compare 
them with models in the literature. Ideal candidates for future searches are 
intrinsically faint stars, such as M dwarfs, at very close distances. 

Subject headings: stars: late-type, extrasolar planetary systems, planetary: at- 
mosphere, radiation mechanism: non-thermal, stars: individual: v And, HD 
46375, 55 Cnc, GJ 436, r Boo, HD 217101 

1. Introduction 

Twenty percent of the 180 known extrasolar planets orbit within 0.1 AU of their parent 
star and have masses comparable to Jupiter's. The successful radial velocity technique (e.g.. 
Mayor & Queloz 1995; Marcy & Butler 1996) provides the Keplerian parameters for these 
'hot Jupiters': minimum mass (Mpsini), orbital period (Porfe), semi-major axis (a), and eccen- 
tricity (e). For the very small number of these that transit their star, the transit lightcurves 
supply the planet radius (Rp), orbital inclination (i), and true mass (M^; e.g. Charbonneau 
et al. 2000, 2006), and in the rare case of a bright host star, spectroscopy of the transit event 
can allow one to measure or constrain elemental and molecular abundance in the planet's 
atmosphere. Our knowledge of a hot Jupiter's atmospheric composition is limited to the case 
of HD 209458 b, whose parent star is bright enough (V=7.65) for 'transit spectrophotome- 
try.' This method provided the first detection of an extrasolar planetary atmosphere through 
observations of Na 1 absorption by the planet's troposphere (Charbonneau ct al. 2002). 

Vidal-Madjar et al. (2003, 2004), using transit spectroscopy, reported a detection of 
ultraviolet absorption by H 1, O I and C II in the upper atmosphere of HD 209458 b out 
to several planetary radii. ^ They inferred that the observed 15% Lya absorption is taking 
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place beyond the planet's Roche hmit, and thus that hydrogen was escaping at a rate of 
~10^° g s~^. This observation motivated several models of the upper atmosphere of hot 
Jupiters that experience the influence of the parent star's extreme ultraviolet (EUV) flux on 
temperature and escape rates. Rapid hydrodynamic escape of hydrogen would potentially 
affect the planet's structure, evolution, and shed light on the mass-radius relationship of 
the known transiting planets (Pont et al. 2005). The first such calculation was presented 
by Lammer et al. (2003), who based their model on Watson et al.'s (1981) hydrodynamic 
treatment in which upper atmospheric temperature and mass loss are coupled. Lammer et 
al. pointed out that even though the effective temperatures of hot Jupiters are low enough 
for the planets to be stable against mass loss, UV-heated exosphere temperatures are much 
higher, promoting escape. They calculated a H escape of ~10^^ g s~^, 9 orders of magnitude 
greater than the Jeans escape at the planet's effective temperature (Sasselov 2003). Baraffe 
et al. (2004) applied this work to evolutionary calculations where they unite irradiated planet 
atmospheres with interior structures. They show that for a given orbital distance there is a 
critical mass below which the planet will undergo rapid runaway evaporation, and though 
extremely improbable, HD 209458 b might be in such a fleeting phase. 

Using a more detailed model, Yelle (2004) calculated that the energy-limited atmo- 
spheric escape rate is approximately proportional to the stellar EUV flux and for HD 209458 
b, is 20 times less than the Lammer et al. value. This difference is ascribed to the depen- 
dence of escape rates on the thermal structure of the upper atmosphere and the previously 
neglected balance between heating by EUV radiation and coohng by escape of H and radi- 
ation by molecules in the homopause at infrared wavelengths. Yelle (2004) showed that the 
near IR emission by Hj]" is an important cooling mechanism in the upper atmosphere of hot 
Jupiters just as it is in Jupiter's hot thermosphere. 

In the following section, we present the predictions of H3" emission from a hot Jupiter 
which encouraged this observing program. In addition to being the first detection of Hg" in 
emission outside of the Solar System,^ direct planetary detection would offer a new view into 
the physical properties of hot Jupiter thermospheres, applying needed constraints to models 
of their thermal structure. We present our search for the Q(1,0) transition of (A = 3.9530 
Hm or 2529.72 cm~^) from the upper atmospheres of six close-in extrasolar planets. This 
wavelength is in the L' passband, where the parent star, relative to optical wavelengths. 



Hubble Space Telescope. Because STIS is no longer operational, the observations cannot be presently 
confirmed. 

^Hg" has been observed in absorption in the interstellar medium (e.g. Goto et al. 2002). Brittain & Rettig 
(2002) announced a detection of Hjj" emission from the protoplanetary disk around HD 141569, but more 
detailed observations by Goto et al. (2005) failed to confirm the result. 



-4- 



is dim making the contrast with planetary emission more favorable. We capitalized on 
the time-dependent Doppler shift of the planet, which can be as large as 150 km s~^ (or 
AA = 0.0020/um at 4/um) between two nights. We differenced spectra to remove the stellar 
photospheric signal and telluric absorption lines, isolating any residual Hj]" emission which 
would stand out as a peak and a velocity shifted dip in the final residual. Section 3 outlines 
our observations and data reduction, while the analysis and resulting luminosity limits are 
presented in Section 4. We summarize our results and discuss the requirements of future 
experiments in Section 5. 



2. H3 emission from hot Jupiters 

Jovian has an equatorial column density of ~10^^ m~^ and is ~100 times more 
abundant at auroral latitudes (Lam et al. 1997). It has been successfully used as a diag- 
nostic of the thermal and chemical state of Jupiter's thermosphere and ionosphere (Stallard 
et al. 2002).^ has a strong ro- vibrational spectrum emitting at 3—4 /xm, correspond- 
ing to excitation temperatures of ~1000 and 750 K, and thereby effectively samphng the 
temperature range of Jupiter's upper atmosphere (Stallard 2001). 

Several groups modeled the formation of H3" in giant planet ionospheres (e.g. Kim et 
al. 1992, Achilleos et al. 1998). A chain reaction starting with the ionization of H2 forms 
the molecular ion in the thermosphere. At high latitudes, this is achieved by collisions with 
energetic electrons funneled down along magnetic field lines: 

H2 + e* ^ H++ e- + e-, 

and more globally, through stellar EUV radiation, 

H2 + hjv ^ H^-F e~. 

An exothermic reaction quickly converts the ion to Hj, 

H+ + H2 ^ H+ + H. 

has a minimum lifetime of about 10 s (Achilleos et al. 1998) and is destroyed by disso- 
ciative recombination: 

H+ + e- ^ H+ + H + e- 
or H+ + e- ^ H + H + H. 



^The solar EUV radiation and magnetospheric charged particles that are absorbed by Jupiter's ther- 
mosphere dissociate or ionize the thermospheric molecules creating a coincident ionosphere (Yelle & Miller 
2004). 
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Jupiter emits ~10^^ W in the emission lines of (Lam et al. 1997). At a distance 
of 10 pc, the resulting flux of SxlQ-^^ W m'^ (=2x10-^0 W m'^ jj,m-^ or 0.03 /xJy if 
placed in an unresolved line at spectral resolution of 36,000) would be undetectable with 
today's telescope technologies. A hot Jupiter (a ~ 0.05 AU) is 100 times closer to its parent 
star and experiences 10'^ times the EUV flux and enhanced magnetic interaction with the 
stellar magnetosphere (Shkolnik et al. 2005). Direct detection of a transiting or even non- 
transiting hot Jupiter atmosphere is possible because a significant fraction of this additional 
energy is re-radiated by narrow lines of molecular coUants, including Hg . As predicted by 
Yelle (2004) for hot Jupiters and measured by Stallard et al. (2002) in our own Jupiter (see 
below), H3 emission is the dominant coolant of the planet's thermosphere. A measurement 
of emission would provide needed information about temperature and ion density in the 
thermosphere of a hot Jupiter. 

Miller et al. (2000) used the Jupiter Ionosphere Model (JIM; Achilleos et al. 1998), a 
three-dimensional, fully coupled model of the Jovian thermosphere/ionosphere, to calculate 
that a Jupiter orbiting the Sun at 0.05 AU would have an Hjj' column density (due to 
irradiation alone) at the sub-solar point of 10^^ m^^, or 1000 times that of Jupiter. The 
predicted emission is ~10^^ W, or a few xlO^^^ W m^^ if the planet is 15 pc away. They 
consider this a lower limit because their model does not take into consideration an extended 
atmosphere as observed for HD 209458 b, a deeper atmosphere for planets that are more 
massive than Jupiter, any increased ionization due to the strong stellar wind, or interactions 
between the magnetospheres of the planet and the star. 

Yelle (2004) presented a 1-dimensional aeronomical calculation representing a global- 
average of the atmospheres of hot Jupiters with a — 0.01, 0.05, and 0.1 AU. The model 
shows that the planets' thermospheres are heated beyond 10^ K by the EUV radiation from 
the solar- type star. The lower thermospheres, between 1 and 1.1 Rp, are primarily cooled by 
Hg radiative emission, despite the fact that the calculated Hg densities are approximately 
a factor of 10 lower than in Jupiter. This is because the increased photoionization both 
creates H^, which is quickly converted to H^, and generates the electrons that destroy the 
H^. Cooling rates vary from 6x10^*^ ergs cm"^ s^^ at 0.1 AU to 3x10^® ergs cm^^ s^^ at 
0.01 AU. For the reference case, modeled after HD 209458 b (a = 0.05 AU, Rp = 1.4Rj), 
Yelle predicted an H^ luminosity of 1.0x10^^ W, an order of magnitude less than Miller et 
al.'s prediction. 

Neither model includes the effects of the stellar gravitational force, stellar radiation 
pressure, or a planetary magnetic field. The first two become important only beyond ~3 Rp, 
however a magnetic field would make a direct contribution to the ionosphere. Though we 
do not yet know how ubiquitous magnetic fields are among extrasolar planets, Shkolnik 



-6- 



et al. (2003, 2005) presented evidence suggesting strong fields {B > Bj — 4.3 G) on two 
hot Jupiters, HD 179949 b and v And b. A planetary field would trap ions, limiting ion 
escape, and cause the precipitation of electrons and ions along magnetic field lines, producing 
polar enhancements of Hj|" on the close-in planets similar to those observed on our Jupiter. 
This might increase the emission from a hot Jupiter, potentially by orders of magnitude 
beyond the predictions. 

3. Observations &c data reduction 

Our six stellar targets have a broad range of spectral types (F7— M2.5) and L' magni- 
tudes (2.8—7.8). Each star hosts either a hot Jupiter- or a hot Neptune-mass inner planet 
with an orbital period of <7 days. Given the uncertainties in the predictions, there was 
a possibility that detectable emission might have come from planets farther away from 
their central star. For this reason, we included three stellar systems that have one or more 
additional outer planets. Table 1 lists each planet's minimum mass, orbital period, and 
semi-major axis, while the stellar properties are given in Table 2. 

We secured two nights of L'-band spectra of the six stars on the 3-m NASA Infrared 
Telescope Facihty's (IRTF) Cryogenic Near-IR Facility Spectrograph (GSHELL) in its high- 
resolution echelle mode during 2005 February 15—17 and August 26 & 28. (August 27 was 
lost to poor weather.) A 0.5" x 30" sht gave a spectral resolution of 36,000 with an average 
resolution element of 2.70 detector pixels (8.3 km s~^). Narrow-band circular variable filters 
(GVF) isolate a single order within a spectral range of 1— 5/im. The 256x160 SBRG InSb 
detector (Tokunaga et al. 1990; Greene et al. 1993) intercepts a narrow wavelength range 
equal to 1/400 of the central wavelength, in this case, 3.948—3.958 /xm. 

Wc nodded the star along the slit with a beam pattern of A-B-B-A to ensure accurate sky 
subtraction. The detector's response becomes non-linear beyond 3000 DN (33,000 e~), which 
forced us to keep all exposures no longer than 70—120 seconds. The total integration times 
and resulting S/N for each star are listed in Table 2. The data reduction was performed 
with standard IRAF'^ routines. We first applied a bad pixel mask to correct for ~250 
"noisy" (non-linear) pixels, and eliminated cosmic ray hits with the cosmicray procedure. 
We produced a nightly mean dark-subtracted, normalized 2-D fiat field to flatten each [A— B] 
and [B— A] frame. Movement of the GVF wheel resulted in asymmetric fringing and imperfect 



^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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flat-fielding (Figure 1, J. Rayner, personal communication). This produced inconsistencies 
between the A and B beams and cannot be avoided with CSHELL since one needs to move 
the CVF wheel to acquire a star, and again to take a fiat field exposure. We therefore treated 
each beam separately throughout the data reduction and analysis, combining only the final 
A- and B-beam differenced spectra (residuals, where the fringing effects are removed) for 
increased S/N. 

We extracted a 1-D spectrum from each sky-subtracted, fiat-ficldcd 2-D image and 
calibrated wavelength with a krypton arc lamp spectrum. We took arc exposures before and 
after each stellar target. Due to the very small spectral range, it is rare for an arc line to fall 
within this span. It was necessary to observe arc lines at the same grating position as used 
for H3 line, move the CVF filter to select three different grating orders, and extrapolate to 
the desired wavelengths using the grating law (relating orders, angles, and wavelengths; see 
CSHELL manual). One-dimensional arc spectra were extracted for each of the two beams 
since the curvature in the arc lines across the detector produced a small difference in the 
wavelength solutions. The final precision of the wavelength solution is 2.5xl0~5 /im. 

We were unable to remove the N2O telluric absorption lines due to the poor S/N of our 
hot standard star spectra. This did not pose a problem for most of the stellar targets when 
comparing two nights of data because the atmospheric N2O lines were stable. However, 
we could not apply heliocentric velocity corrections with the telluric lines in the spectra, 
or "shift & add'' the individual exposures to accommodate for potential broadening of the 
Hjj" line due to the planet's orbital motion within the nightly span of observations. These 
velocities are for the most part negligible at this spectral resolution. The heliocentric (and 
bary centric) velocities differ from night to night by less than 3 km s~^. The broadening of 
any planetary emission due to orbital motion in a couple of hours is usually <8 km s~^, or one 
resolution element. Broadening due to the planet's rotation is also insignificant since these 
hot Jupiters are hkely tidally locked to their stars, resulting in rotation rates of <1 km s~^. 

In the end, spectra of each target star on a given night were combined to produce 
the final A- and B-beam spectra. The sought-after H3 emission would be greatly Doppler 
shifted (150 km s~^) and stand out as a peak and a velocity shifted dip in the mean residual 

{ — ^[{■^Nightl — ^Night2) + (Bmghtl — ^mght2)])- 

4. Results 

L'-band spectra of v And, one of our brightest targets, are displayed in Figure 2. We plot 
the normahzed A- and B-beam spectra taken on 2005 August 26 and 28 in the top two panels 
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with their corresponding residuals in the third panel. A Fourier analysis revealed a peak 
frequency of 575 /xm"-*^ (0.00244 cycles pixel"^) in both the v And spectra and their residuals. 
We removed this frequency from the residuals as it is clearly due to small variability in telluric 
absorption. The final mean residual is shown in the fourth panel of Figure 2. The boxed 
region spans the wavelength range where the Doppler-shifted emission would appear 
if detected. The RMS within this region is 0.38% of the stellar photosphere signal. The 
residual RMS is converted into an Hg luminosity limit using the star's L' magnitude and 
distance. We derived the L' magnitudes from the K magnitudes of the Two Micron All 
Sky Survey (2MASS, Skrutskie et al. 2006) and intrinsic colors from Bessell & Brett (1988). 
For V And, the RMS translates into a luminosity limit of 1.3x10^^ W. The RMS of the 
intranight residuals for the v And spectra is 0.2%. This demonstrates the level of stability 
of the spectrograph on short timescales as well as the reliability of the data reduction and 
analysis. 

The normalized spectra of r Boo and 55 Cnc are presented in Figures 3 and 4. For 
these lower S/N spectra the RMS values are 0.92% and 0.98%, respectively, corresponding to 
luminosity limits of 2.2x10^^ and 1.2x10^^ W. In Figure 5, we plot the RMS of the night-to- 
night residuals for the six targets against the theoretical noise limit, (S/N)~^. Observations 
of V And reach this limit but as the S/N decreases, the RMS deviates from the line. This is 
mainly due to systematic noise in the CSHELL detector, including the "noisy" pixels, which 
are more difficult to identify and remove with lower S/N data. The upper limits for all the 
targets are listed in Table 2 and plotted against absolute L' magnitude in Figure 6. Clearly, 
higher S/N data is not the sole contributor to a more stringent luminosity hmit. The 
lowest limit is set by the faintest and lowest S/N star because it is the closest target in our 
sample, GJ 436 (SpT = M2.5). GJ 436, at a distance of only 10.2 pc from Earth, has an 
emission limit of 6.3x10^^ W, thus approaching Miller et al.'s (2000) estimate of -10^^ W. 



5. Summary and discussion 

We searched for the Q(1,0) transition of the H3 molecule at 3.953 //m from the thermo- 

spheres of close-in giant planets around six stars and exploited the Doppler shift due to the 
planet's orbital motion to subtract the stellar photosphere and telluric lines. The majority 
of our upper limits on luminosity are higher than predictions by Miller et al. (2000) and 
Yelle (2004). GJ 436, our faintest, nearest and lowcst-S/N target, set the most stringent 
limit of 6.3x10^'' W, comparable to Miller et al."s most conservative estimate, but far higher 
than that of Yelle (2004). Our limits suggest that non-radiation effects (e.g. magnetospheric 
heating) do not dramatically enhance H3 emission from these planets. 
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The brighter (but more distant) stars in our survey offer the least promise for improve- 
ment. The hmit on the emission from the planets orbiting v And, our brightest target 
with the highest S/N spectra, is 1.3x10^^ W. To achieve Miller ct al.'s hmit of ^lO^'^ W, 
a S/N of 3400 pixel"^ is necessary, and 10 times that amount to reach Yelle's limit. This 
would require a prohibitively long exposure time considering that we reached a S/N of only 
300 in 2 hours of integration. We compare these requirements with the capabilities of NIR- 
SPEC, the near-IR spectrograph mounted on the Keck II 10-m telescope: A S/N of 3400 
pixel"^ at NIRSPEC's comparable spectral resolution requires 7 hours of integration time, 
but systematic effects would inhibit achieving such a high S/N. 

The situation is better for less luminous but closer stars for which the contrast ratio 
between planetary H3 emission and the stellar photosphere is higher. The only such star 
in our sample, GJ 436, requires a relatively low S/N of ~100 pixel"^ to achieve a limit of 
10^^ W. This can be reached in 10 hours with CSHELL and in only one hour with NIRSPEC. 
Yelle's limit of ~10^^ W remains unfeasible even with Keck, the largest optical telescope in 
the world. 

Yelle (2004) suggested that 'occultation spectroscopy', as discussed by Richardson et 
al. (2003), may be the best way to observe H3 in planetary atmospheres. HD 209458 b 
(Charbonneau et al. 2000, Henry et al. 2000), TrES-1 (Alonso et al. 2004) and HD 189733 b 
(Bouchy et al. 2005) are currently the only known planets for which this is practical. Given 
that our intranight residual RMS is half that of the RMS between nights, comparing spectra 
before, during and after secondary eclipse is an attractive option. However, one would not 
get the required S/N for these stars during the ~2 hours of a single eclipse event, meaning 
cumulative data over several transits would be necessary. A Jupiter-sized planet transiting 
a nearby M dwarf would be the most promising target for a future H3 search, though such 
a planet has yet to be discovered. 

We thank John Rayner, Alan Tokunaga, Tim Brown, and Sean Brittain for helpful 
discussions. E.S. acknowledges the financial support of the NRG's Research Associateship 
Program. E.G. acknowledges support by the NASA Terrestrial Planet Finder Foundation 

Science Program. This research has made use of the VizicR catalogue access tool, GDS, 
Strasbourg, France, and the Extrasolar Planet Catalog maintained by Jean Schneider. 
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Table 1. Planetary systems observed 



Star 


Planet 


Msini 


^ orb 


a 


Reference 






(Mj) 


(days) 


(AU) 




GJ 436 


b 


0.07 


2.64 


0.03 


Butler ct al. 2004 


55 One 


e 


0.045 


2.81 


0.038 


McArthur et al. 2004 




b 


0.784 


14.67 


0.115 


Butler et al. 1997 




c 


0.217 


43.93 


0.24 


Marcy et al. 2002 




d 


3.92 


4517.4 


5.257 


Marcy et al. 2002 


V And 


b 


0.69 


4.617 


0.059 


Butler et al. 1997 




c 


1.89 


241.5 


0.829 


Butler et al. 1999 




d 


3.75 


1284 


2.53 


Butler et al. 1999 


T Boo 


b 


4.13 


3.3135 


0.046 


Butler et al. 1997 


HD 46375 


b 


0.249 


3.024 


0.041 


Marcy et al. 2000 


HD 217107 


b 


1.37 


7.1269 


0.074 


Fischer ct al. 1999 




c 


2.1 


3150 


4.3 


Vogt et al. 2000 
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Table 2. Targets and observations 



Star 


SpT 


Dist. 
(pc) 






HJD - Night 1 
HJD - Night 2 


Total Exp. 
Time (s) 


Total S/N 
(pixel^^) 


RMS of 
Residual 


Emission 
Limit (W) 


GJ 436 


M2.5 


10.2 


6.07 


5.80 


2453417.97 


6480 


30 


0.0493 


6.3E+17 












2453418.98 


5760 


40 






55 Cnc 


G8 V 


13.4 


4.02 


3.96 


2453417.87 


6360 


310 


0.0098 


1.2E+18 












2453419.85 


960 


100 






V And 


F8 V 


13.47 


2.86 


2.82 


2453610.10 


7280 


310 


0.0038 


1.3E+18 












2453612.06 


7520 


370 






T Boo 


F7 V 


15 


3.51 


3.47 


2453418.06 


3840 


270 


0.0092 


2.2E+18 












2453420.05 


2520 


180 






HD 46375 


Kl IV 


33.4 


7.85 


7.79 


2453417.81 


2640 


30 


0.0462 


l.OE+18 












2453419.75 


7440 


30 






HD 217107 


G8 IV 


37 


4.54 


4.49 


2453609.94 


4320 


160 


0.0196 


l.lE+19 












2453611.91 


4080 


80 







''K magnitudes taken from the Two Micron All Sky Survey (2MASS) via VizieR. 
''Calculated from K magnitudes and intrinsic colors from Bessell & Brett (1988). 
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Fig. 1. — Image of the normalized flat field. The lower and upper strips represent the 
locations of the A and B beams, respectively. The fringing that induces variations in the 
two beams is evident. 
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Fig. 2. — L' spectra of v And centered on tlie 3.953-/im line. Two nights of CSHELL's 
A-beam and B-beam are plotted in the first two panels. The y-axes are normalized intensity. 
The third panel shows the residuals of the two beams. The periodic structure in the beams 
reflects the change in the N2O telluric lines and is removed in the bottom panel which plots 
the mean residual. The boxed region spans the wavelength range where the Doppler-shifted 
H3 line would be and has an RMS of 0.38%. 
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Fig. 3. — L' spectra of r Boo centered on tlie 3.953-/im line. Two nights of CSHELL's 
A-beam and B-beam are plotted in the first two panels. The y-axes are normalized intensity. 
The third panel shows the mean residual of the two beams as discussed in the text. The 
boxed region spans the wavelength range where the Doppler-shifted line would be and 
has an RMS of 0.92%. 
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Fig. 4. — Spectra and tlie mean residual of 55 Cnc. See caption to Figure 3. 
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Fig. 5. — Residual RMS plotted against (S/N) ^ for each target star. The red line represents 
the theortical S/N limit, v And is the only star for which we reach this limit. 
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Residual RMS converted into equivalent H3 luminosity limit plotted against ab- 
of the star. 



